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Abstract 

Acyclic networks are a class of complex networks in which links are directed and don't have closed 
loops. Here we present an algorithm for transforming an ordinary undirected complex network into 
an acyclic one. Further analysis of an acyclic network allows finding structural properties of the 
network. With our approach one can find the communities and key nodes in complex networks. 
Also we propose a new parameter of complex networks which can mark most vulnerable nodes of 
the system. The proposed algorithm can be applied to finding communities and bottlenecks in 
general complex networks. 
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I. INTRODUCTION 



Acyclic networks are a class of complex networks in which all the links are directed and 
don't have closed loops. As an example the network of food chains is the acyclic one: 
edges are directed from anode representing a predator to the node representing its prey 
Another example of an acyclic complex network is the network of citation: this is a directed 
network in which every edge has a direction {2]. 



Basic research in the field of acyclic complex network is focused on networks with a hi- 
erarchical structure. Systems, organization of which displays causal asymmetry constraints 
can be described in terms of directed links using a discrete set of arbitrary units, the evo- 
lutionary history of a given collection of events or the chart of computational states visited 
along a complex simulation. Such a set of paths and links defines the acyclic network 3|-[5]. 
Acyclic complex networks are widely used for different physical problems: from protein fold- 
ing to sociophysics 6, 3]. One of the fundamental problems of the complex networks is how 
to divide network into its constituent communities. In the complex networks communities 
are the groups of nodes that are densely connected amongst themselves while being sparsely 
connected to the rest of the network. There are different methods to detect communities 
based on calculation of betweenness, loops of vertexes, analysis of adjacency matrices of the 
network js-16]. 



In this paper we propose a method to map a general complex network into acyclic graph, 
arranging them in strong hierarchical order. Mapping a complex networks into the acyclic 
one opens new opportunities to analyze them which were not present in the usual complex 
networks. The calculation of these new parameters allows more detailed analysis of the 
structural characteristics of a complex network. With the proposed transformation of the 
network it's possible to detect most vulnerable nodes of the structure and detect cluster 
composition of the network. With our approach one can find the communities, hubs and 
key nodes in complex networks. Also we propose a new parameter of complex networks 
which can mark most vulnerable nodes of the system. The proposed algorithm gives a clear 
picture of an investigated systems, describing some structure details of the complex network. 
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II. METHODS 



An acyclic network is a graph in which the links have direction and do not form loops. 
In order to transform a complex network into acyclic one, it is proposed to apply to two 
nodes of the network the potential difference that will flow from its current. The direction 
of current flow through the links will set the direction of the links, so in this way the edges 
don't form a loop. If the network nodes are not connected with the others, then no current 
will flow through them and acyclicity will not be broken. 

To apply current we choose two nodes which we call input and output and put their 
Voltage to 1 and respectively. The resistance of each link is inverse to its weight: 

Rij = ±- (l) 

Then the Kirchhoff equation can be written as: 

Vinput = 1 (2) 
y output = (3) 
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where hi is node degree and a^- is the adjacency matrix of the network. In general solving 
the Kirchhoff equations takes 0(n 3 ), but there was proposed the method when the voltage 
values can be obtained in linear time Q. After calculating all the potentials we arranged 
all the nodes in the network in the acyclic sequence according to their potential (See Fig.l 
for example). To characterize the whole network we chose two nodes, which were the most 
important for our study. As an example they can be two servers of computer network, if 
we are interested in the channel capacity between them, or^they can represent folded and 
unfolded states of protein in conformation space network 



After applying the potential difference we obtained the sequence of nodes and we could 
relabel them according to their voltage. Once the nodes are ranked, we can work with 



transformed network. In jo] there was introduced the notion of setting the node flux Aj and 
the flux between two selected nodes /Xj, which is defined as: 



i-l 



3=1 3=1 
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where /c m and & OMi stand for the links going out from the node and going in respectively 
and kj is the index of the node in acyclic sequense i.e. node ranking according to its potential. 
Parameter Aj characterizes the amount of links that connect the nodes to the i-th and after 
it. (See Fig. 2 for visual explanation) 



III. RESULTS 

For analysis we created Watts-Strogatz and Erdos-Renyi networks with number of nodes 
from 200 to 500 and the following probabilities of connection 1% for Watts-Strogatz and 



1%,2%,3% for Erdos-Renyi networks 
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To facilitate calculation in our model, all 
nonzero edges have resistance R = 1. For a single run (potential difference was applied 
only to one pair of nodes) we observe different behavior for Watts-Strogatz and Erdos-Renyi 
(Fig.3). In Watts-Strogatz networks the number of links is smaller comparing to Erdos-Renyi 
and additional links which were randomly put in the network create small separated clusters 
that can be observed. This effect becomes stronger with the increasing of the probability of 
adding random links, which is crucial for creating subclusters in this type of networks. 

For Erdos-Renyi type of networks a single-well structure was observed and we can con- 
clude that the connectivity is homogeneous and there are no well defined clusters. The 
analysis show no principal difference for all possible pairs in Erdos-Renyi networks. On the 
other hand, small-world networks' node flux demonstrates splitting the network into many 
small clusters. So we can make a conclusion that with the observed node flux behavior it's 
possible to detect different clusters in the network. Although, if both key nodes are taken 
from the same cluster one will get single-well picture, so single-run calculation is not enough 
to describe the whole structure of the system. 

For the further observation a network consisting of two Erdos-Renyi subnetworks was 
created, which connected a small number of links (comparing to the quantity of links in the 
subnetworks). Calculated node flux (Fig. 4) accurately describes the structure of the created 
system and size of its components. We analyzed this network to check how sensible was the 
obtained structure with different pairs of input and output nodes. For the first calculation 



the nodes belonging to the different subnetworks were taken as input and output and then we 
took pairs of nodes by random. In both cases we took 1000 pairs of nodes. Averaged results 
for both measurements are shown in Fig.5. For pairs of nodes from different subnetworks it is 
easy to observe the structure of the network. For the other case picture is blurred, because 
if both input and output were taken from the same cluster they don't give the picture 
of the real structure. This problem, however, can be solved by using sufficient statistical 
data. Of course, the precision and accuracy of the obtained structure in case of many pairs 
measurement are strongly dependent on particular network. Thus, this approach allows 
detecting the cluster structure of the network and determine belonging of certain nodes to 
cluster in the network. 

It is clear that in case of many communities the general picture of Aj will not show the 
real structure of the network. E.g. if there are two "parallel communities" with the respect 
to the input community, their nodes will be mixed. But in this case our approach can define 
the community of input node with the same quality as in the case of two communities. So 
one should perform the algorithm for the i node, get the list Q of the nodes belonging to 
the same community and than run algorithm with another input not belonging to the C{. 
To test our approach we performed the set of benchmarks proposed by Fortunato et al. 



2l|. One of the main parameters in this benchmark is the mixing parameter,//, which is 
defined as fraction of the links which node shares with another communities(i.e. if yU=0.8 it 
means that 80% of the links of the node go inside its community and 20% are shared with 
other communities). Degree distribution of the node and community distribution are taken 
from the power law distribution with exponents 7 and (3 respectively .It's easy to see that 
/i = 0.5 or higher our approach will fail because there is no boundary between community 
of the node and the rest of the network. We tested our approach on the networks with the 
number of nodes N=1000,/3 = 1,2 and 7 = 2, 3(Fig. 7). The algorithm shows good results for 
values of < 0.3 but for further values its precision decreases dramatically. The results are 
slightly clearer if the average degree of the node is increased because the boundaries between 
the communities become more pronounced. The comparison of most popular algorithms was 



done by 



0, E 



m- 



jancichinetti and Fortunato [22j . Although the results of some popular algorithms 
30] are better, the method proposed here can be performed in linear time which is 
its main advantage. 

We also introduced the vulnerability of the single node A i; that can be used to identify 



the most vulnerable nodes of the network. Parameter characterizes the ratio of flux which 
is flowing through the given node (blue edges) and its total flux (blue and red edges). Aj is 
defined as the ratio of links that go across the given node and its total degree: 



^ _ fa — \ + Hi+l — \ _ ki ^ 

In these terms the proposed parameter can describe the importance of the chosen node 
by indicating the value of the flux which is going from one part of the system to another 
precisely through the given node. Apparently if the node is a bottleneck between two or 
more clusters its vulnerability will be maximum, because all the flow between two parts of 
the system will be flowing only through the link which is connecting two boundary nodes. 
To test our hypothesis we calculated the vulnerability for two different networks: normal 
Erdos-Renyi network and network which has two Erdos-Renyi networks connected by small 
number of links. In the first case Aj is maximum at the beginning and at the end of the 
acyclic sequence, which can be explained by the small number of links connected to the 
input and output nodes; in the case of two clusters we observed the same regions of high Aj 
values but also another maximum region in the place of connection of two clusters, which 
means that the nodes in this region are highly loaded and among important for connection 
of input and output nodes. 



IV. CONCLUSION 

In this paper we present a method to transform general complex network into the acyclic 
one. The proposed method can be used to determine cluster structure of the network and the 
size of its components. According to our measurements the node-flux shows different trend 
for Erdos-Renyi and Watts-Strogatz types of complex network in single-run calculations, so 
we can conclude that it's possible to determine the type of the network. We performed the 
set of the benchmarks which shows that our method can be used to detect communities in 
linear time. The proposed parameter vulnerability can be used to detect key-nodes of the 
network, especially bottleneck-type networks, not only so-called hubs, which can be detected 
by measuring average node flux and vulnerability. 
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FIG. 1. Scheme of the transforming of a complex network into an acyclic one. The link 034 was 
removed because the voltages of the nodes are equal. 




FIG. 2. Visual representation of the node flux, and the flux between the nodes of a complex 
network. 



VI. FIGURES 




FIG. 3. Typical plots of Watts-Strogatz (a) and Erdos-Renyi (b) networks flux, where i is index 
of node in the acyclic sequence. 
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FIG. 4. Flux of the network, which has two Erdos-Renyi subnetworks connected by small number of 
links. Each subnetwork is represented by bell-like flux which is typical for Erdos-Renyi networks. i 
is index of node in the acyclic sequence. 
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FIG. 5. Cut with input and output in two different clusters (a) and with randomly chosen input 
and output (b). 
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FIG. 6. Vulnerability characterizes the ratio of flux which is flowing through the given node (blue 
edges) and total flux (blue and red edges). 
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FIG. 7. Results of the benchmark of the random networks with N=1000. The data are presented 
when average node degree, < k > is equal 10,15,25. 
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FIG. 8. In the figure (a) is shown the vulnerability of a network which has only one Erdos-Renyi 
cluster, in the figure (b) the network consists of two Erdos-Renyi clusters with a small bridge 
between them. % is index of node in the acyclic sequence. 
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